Measurement of CNGS muon neutrino speed with Borexino 
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Abstract 

We have measured the speed of muon neutrinos with the Borexino detector using short-bunch CNGS beams. The 
final result for the difference in time-of-flight between a <E>=17 GeV muon neutrino and a particle moving at the 
speed of light in vacuum is 5t - 0.8 + 0.7 stat + 2.9 sys ns, well consistent with zero. 

Keywords: neutrino speed, special relativity, GPS time-link 
PACS: 03.30.+p, 14.60.Lm 



Introduction 

This paper describes a precise measurement of the 
speed of CNGS H] muon neutrinos made with the 
Borexino [0 [3] H) detector at the Laboratori Nazionali 
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5 del Gran Sasso (LNGS) in Italy. 

e CNGS neutrinos travel about 730 km in matter with 

7 one of the highest relativistic y factors ever artificially 

s produced. The neutrino mass is at most 2 eV/c 2 or 
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9 possibly much less, while the CNGS average beam en- 

10 ergy is 17 GeV, so y is always >10 n , much bigger than 

11 that obtained in any charged particle beam. A test of 

12 Special Relativity with these particles is therefore mean- 

13 ingful. Besides, the measurement may also put an upper 

14 limit on the effect of non-standard propagation of neu- 

15 trinos in matter. 

16 This effort was also motivated by the claim, made by 

17 the Opera Collaboration in Sep. 2011 [5], that CNGS 

18 neutrinos travel faster than the speed of light in vacuum. 

19 This claim, however, was later withdrawn |6|. 

20 If the mass of the heaviest neutrino is assumed to be 2 

21 eV/c 2 , the best direct limit on a neutrino mass, then the 

22 relativistic velocity of a 17 GeV neutrino should satisfy 

23 |v - c|/c < 1CT 19 . Cosmological measurements Q and 

24 neutrinoless double beta decay measurements |8| give a 

25 mass limit an order of magnitude smaller, implying an 

26 even smaller constraint. However, theories with extra 

27 dimensions [ 9 1 predict apparent velocities different than 

28 the speed of light. Some of these theories IfTOl [TT1 [T2l 

29 allow |v-c|/c ^ 10~ 4 at neutrino energies of a few GeV. 

30 The present work has the sensitivity to test these theo- 

31 ries and put limits on them. 

32 The MINOS collaboration has already performed a 

33 measurement with the NuMI beam, yielding a result 

34 compatible with c with an uncertainty of ^10~ 4 IT51 . 

35 A similar measurement was also recently done by the 

36 Icarus collaboration [ 14 1 with a small sample of events 

37 collected in the Nov. 201 1 CNGS run. 

38 We have been collecting CNGS events since May 

39 2007, when the Borexino experiment started taking 

40 data. However, the triggering and time tagging systems 

41 originally designed for Borexino 021, while perfectly 

42 adequate for solar neutrino physics lfl6l [T71 [THl and su- 

43 pernovae detection, were not sufficiently precise to al- 

44 low an interesting measurement of the neutrino velocity. 

45 We have therefore designed and installed a new facility 

46 capable of achieving an uncertainty of the order of a few 

47 ns. This new facility is based on a small-jitter analogue 

48 trigger, a geodetic GPS receiver, and a GPS-disciplined 

49 atomic Rb clock. The reader interested in the techni- 

50 cal details and the calibration of this system may refer 

51 to fT9l. 

52 The Borexino detector is a high-purity liquid scin- 

53 tillator calorimeter (within a Stainless Steel Sphere of 

54 a 1300 m 3 , SSS) shielded by a large Water Tank (WT, 

55 -3500 m 3 ) that serves also as a muon detector. It is in- 

56 stalled in Hall C of the LNGS at a depth of 3800 m.w.e. 

57 CNGS muon neutrinos are detected via charged cur- 

58 rent interactions that mostly occur in the rock upstream 

59 the detector. Internal events exist, but they cannot be 
eo easily disentangled from the crossing muons. Their 



61 number, however, is small 5%). 

62 The kinematics of the muons does not affect the pre- 

63 cision of the measurement. The difference in the time- 

64 of-flight of a « 10 GeV muon with respect to a neutrino 
es traveling a length of about 50 m (the approximate aver- 
ee age distance traveled by such a muon in the rock plus 
67 the distance of the Borexino detector from the North 
ea wall of the Hall C) is less than 0. 1 ns and can therefore 

69 be neglected. The same argument applies to the pions 

70 generated in the CNGS target at CERN. 

71 Muons can be detected by Borexino using the Wa- 

72 ter Tank (Cherenkov light) or the scintillator detector. 

73 In this analysis we have used only events detected by 

74 the latter since its precision is higher. Although only 

75 the core 270 t of PC-PPO scintillator are normally used 

76 for solar and geo-neutrino physics |20|, the sensitive 

77 mass for CNGS muon neutrinos is made by the whole 

78 a; 1300 m 3 of liquid scintillator and buffer liquid be- 

79 cause, though quenched by a small amount of DMP (3 
so g/L), the buffer liquid light yield is sufficient to detect a 
si muon (^2 MeV/cm, equivalent in the buffer to approxi- 

82 mately ^50 p.e./cm) . 

83 The cross-section of the WT is 266 m 2 while that of 

84 the SSS is 147 m 2 . These large areas yield a large num- 

85 ber of events: in May 2012 we have collected a total of 
se 291 events, 144 crossing the SSS and 147 crossing the 
87 WT only. 

ss The distance of the SSS from the CNGS proton tar- 
so get was determined by performing a complete geodetic 

90 determination of the position in the IGS08 Reference 

91 Frame of the underground laboratory and of the Borex- 

92 ino detector in particular. The measurement was done 

93 in collaboration with LVD and ICARUS. The position 

94 of the CERN target was not measured again. We rely on 

95 the precise knowledge of the location of each element 

96 of the accelerator already available at CERN. 

97 The paper is structured as follows: in next section we 

98 describe the measurement of the distance flown by the 

99 neutrinos; we then describe the analysis of the data col- 

100 lected in May 2012; for completeness, we also briefly 

101 report about the result obtained (with much less preci- 

102 sion) in Oct. -Nov 201 1. Finally, we report our result. 

103 Geodesy 

104 The distance between the Borexino reference point at 

105 LNGS and the target at CERN has been computed based 

106 on a geodetic survey at Gran Sasso Laboratories and the 

107 existing target co-ordinates at CERN. 

ioa The geodetic campaign at LNGS was carried out in 

109 May 2012 and was performed in two steps. In the first 
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Figure 1: Top: Cumulative distribution of pots delivered by CNGS 153 

and the equivalent number of pots collected during May 2012 run. 154 

The corresponding live-time is 97%. Bottom: The pot delivered as a 155 
function of time during May 2012 run. 



step two GPS networks, one local and one regional, 
were established. The regional network used in this 
work consists of 32 GPS permanent stations located be- 
tween CERN and Gran Sasso (including two antennas at 
CERN and one at LNGS) whose position in the IGS08 
reference frame has been precisely estimated by adjust- 
ing (fitting) two weeks of GPS data (Bernese software 
ED, GPS weeks 1684 and 1685). 

A local 9 points GPS network has been monumented 
(installed) in the Gran Sasso area and framed to the 
regional one through 3 common stations using a four 
days-24 hours campaign. In the second step, based on 
the previously estimated GPS points, a high precision 
traverse has been measured along the Gran Sasso tun- 
nel highway (10.5 km). This allowed the co-ordinates 
estimation of the Borexino reference point inside the 
LNGS. 

The IGS08 co-ordinates of this point and the re- 
lated standard deviations are: (X= 4582157.919 m, 
Y= 1106469.341 m, Z=4283636.931 m); (o^ =0.020 
m, <x y =0.021 m, cr z =0.015 m). The co-ordinates and 
the precision of the target point at CERN, coded as 
TT41_T_40S, have been supplied by the CERN geode- 
tic team. After datum shift to IGS08, assuming tar- 
get point co-ordinates precision of 0.030 m, the geo- 
metric distance between the Borexino reference point at 
LNGS and the target at CERN has been estimated in 
730472.082 ± 0.038 m. The error should be considered 
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as 1 cr. 

Analysis of May 2012 bunched beam data 

A special beam was set up at CERN in May 2012, 
optimized for the neutrino speed measurement. 

The main parameters of this beam are: narrow 
bunches (cr ^2 ns), 16 bunches per batch with a bunch 
separation of ^100 ns, 4 batches per extraction sepa- 
rated by ^300 ns; the bunch intensity is ^10 n protons; 
one extraction per CNGS cycle, with a cycle length of 
13.2 s. 

We have taken data from May 10''' until May 24"\ 
2012. Fig. [T]top shows the cumulative distribution of 
the delivered protons on target (pots) and the equivalent 
pots collected by Borexino. The total pots delivered was 
19.44 10 16 , out of which 18.88 10 16 were recorded. The 
data taking efficiency during this period was 97%. The 
beam delivery was very stable (see Fig. [T]bottom). The 
missing data on May 23 is due to SPS maintenance. 

As already mentioned in the introduction, we have 
collected a total of 291 events of which 144 crossed the 
SSS. The events were selected by requiring them to be 
on-time within a window of 100 fis with respect to the 
closest CNGS GPS tag plus 2.439 ms, the nominal neu- 
trino time-of-flight. 

This large statistics allows us to apply stringent qual- 
ity cuts and select the best data sample for the measure- 
ment. We have chosen not to use the WT only events, 
because the large size of the tank and the very complex 
internal geometry makes the full simulation of the light 
propagation and collection hard. The event sample was 
further reduced by the requirements of the High Pre- 
cision Timing Facility (HPTF, 1(191 ). which is a purely 
analogue trigger with a threshold of 100 mV, which cor- 
responds to s;800 p.e., i.e. a;8 cm of muon track. This 
threshold reduces the acceptance of muons when they 
cross a small segment of the SSS. 

Borexino has developed muon reconstruction soft- 
ware capable of determining the location of the entrance 
point of a muon in the SSS with a precision of about 50 
cm. This makes it possible to correct for the spherical 
shape of the detector which of course implies a different 
time-of-flight for different entrance points. 

This correction reduces the size of the data sample 
because some events cannot be properly reconstructed, 
but narrows the time distribution significantly, improv- 
ing the quality of the measurement. After this final re- 
duction, our data set consists of 62 CNGS events, more 
than enough for a precision measurement. 

The HPTF provides the GPS time of each trigger. 
The time-link between the CERN tagging system and 
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its calibration is available through public CERN reports 



Figure 2: Time distribution of the 98 inner detector events before fold- 
ing of the 64 peaks and before the muon cuts. The horizontal axis is 
the time distance with respect to the first sample of the waveform dig- 
itizer described in the text. The distribution is flat and its total width 
is consistent with the expected 6900 ns. 



the Borexino HPTF at Gran Sasso has been computed 
and calibrated by collaborators of the Italian and Span- 
ish institutes of Metrology (INRIM and ROA, respec- 
tively) using Precise Point Positioning (PPP) El and 
P3 J23 1 algorithms, with latter implementing an All-In- 
View approach [24J. The calibration of the time link 
has been achieved with an uncertainty of 1.1 ns, being 
this the sum of a systematic uncertainty inherent to GPS 
system and available algorithms with an additional sta- 
tistical uncertainty of 0.3 ns and 1 ns for the PPP and 
P3 algorithms, respectively. The procedure adopted to 
calibrate the HPTF itself and the timelink to CERN is 
fully described in fT31 . 

Fig. [2] shows the time distribution of the collected 
events before the folding of the peaks described below. 
The 98 events in the plot are those that survive a charge 
cut (>800 p.e., corresponding to approximately 2 MeV 
equivalent energy deposited in the scintillator) and that 
have valid data in the CNGS database. 

The time structure and the intensity of the proton 
bunches are recorded shot by shot through a Beam Cur- 
rent Transformer (BCTfllocated before the CNGS neu- 
trino target (743.391 + 0.002 m upstream). 

The signal is acquired by a fast (1 Gs/s) digitizer and 
stored in a database. We call this "proton waveform." 
Every proton waveform is time-tagged with respect to 
the CERN SPS timing system. The GPS time stamp 
marks the time of the first sample of the digitized signal. 
A detailed description of the CERN timing system and 
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The time stamp associated with each waveform, after 
correction by some delays, defines the start time of the 
velocity measurement. Some of these delays are fixed 
in time while some of them are time dependent. 

The arrival time in Borexino is provided by the HPTF 
and the time-of-flight is obtained by the difference of 
the two time tags plus a set of corrections due to: 
CERN-HPTF time link computed with PPP/P3 algo- 
rithms, fixed delays at CERN and Gran Sasso, and light 
propagation and trigger formation in Borexino (event 
dependent). 

The event distribution in Fig. [2]is folded into a single 
one by using the information provided by the BCT dig- 
itized waveforms. After the corrections due to known 
instrumental delays, the geometrical correction that de- 
pends on the muon entrance point, the subtraction of the 
expected time-of-flight at speed c, and a 2.2 ns (posi- 
tive) correction due to the Sagnac effect we obtain the 
distribution shown in Fig. [3] The distribution of the fi- 
nal sample of 62 events is well centered around zero and 
has a width of 4.9 ns. 

The width of the time-of-flight distribution is well un- 
derstood. The light propagation in the detector, the col- 
lection of the light by the PMTs, and the formation of 
the trigger occur at a time that depends on the impact 
parameter of the muon track. Although we correct for 
the muon entrance point, the resolution of the recon- 
struction and the fluctuations in the propagation and in 
the collection of the light yield a distribution of finite 
width. We have therefore carefully simulated the broad- 
ening of the distribution due to the known effects with 
both the Geant-4 simulation and with a dedicated fast 
Monte Carlo. The result of both simulations is in good 
agreement with the data. The width of the arrival time 
distribution of the events from Monte Carlo, including 
light propagation, trigger formation, the jitter of some 
electronics modules, and the intrinsic 2 ns width of the 
bunches is shown in Fig. [3] overlaid with data. This 
Monte Carlo predicts a width of 5.1 ns, in very good 
agreement with the observed width of 4.9 ns. 
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259 Analysis of Oct.-Nov. 2011 bunched beam data 

260 In the period from Oct. 21" until Nov. 6 th , 201 1 the 

261 standard CNGS beam was modified to produce 4 short 

262 bunches (cr ~2 ns) separated by 524 ns. The intensity 

263 was 2.5 10 1 1 protons per bunch. In 16 days of operations 

264 CERN delivered 7.43 10 16 pots which we collected with 

26 5 an overall efficiency of 93% (6.87 10 16 pots). 
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Figure 3: Final distribution of the difference between the neutrino 
time-of-flight and that expected for a particle moving at speed c (data 
points). The mean value is consistent with zero and the width agrees 
with Monte Carlo simulation of known time jitters (yellow filled his- 
togram). The small asymmetry is predicted by Monte Carlo and is due 
to the spherical shape of the Borexino detector. The asymmetry does 
not affect the estimation of the central value by more than 200 ps. The 
mean value and its statistical error are obtained from the gaussian fit. 
The systematic error is explained in the text. D is the distance between 
the proton target at CERN and the center of the Borexino detector. 
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Figure 4: Top: the time distribution of the CNGS neutrino candidates 
selected with the Inner Detector only from Oct. -Nov. 2011 bunched 
beam data. The four peaks separated by 524 ns are clearly visible. 
Bottom: the time-of-flight difference with respect to the speed of 
light after all corrections have been applied, including the Sagnac ef- 
fect. The result is consistent with zero. 



266 We took data with the standard triggering and GPS 

267 tagging systems, so the precision of this measurement 

268 is much lower than the one obtained in May 2012. In 

269 particular, the standard Borexino trigger has a relatively 

270 large intrinsic jitter of 32 ns due to the combined ef- 

271 feet of a set of FPGAs and a DSP which limit the final 

272 precision. Besides, the absolute calibration of the GPS 

273 receiver delay in Hall C (ESAT-100 slave) has an uncer- 

274 tainty of a few ns. 

275 The initial sample was made of 116 on-time events 

276 selected using both detectors in a time window of 16 

277 fis around the nominal time-of-flight value. After se- 

278 lection of internal detector events, pile-up removal and 

279 quality cuts a sample of 36 events was used to measure 

280 the speed of neutrinos. 

281 We have applied to the final sample an analysis pro- 

282 cedure very similar to that performed for the May 2012. 

283 The result is shown in Fig. [4] As it is clear from the 

284 figure, the central value is well compatible with zero, 

285 but the error is much larger than the one obtained with 

286 May 2012 data. The statistical error comes from the 

287 width of the distribution. The systematic error, besides 

288 the contributions described in the next section, is bigger 

289 because of the poor knowledge of the ESAT-100 GPS 

290 receiver calibration and internal delays. 



291 Results 

292 The final result obtained in May 2012 for the time- 

293 of-flight difference of <E>=17 GeV muon neutrinos 

294 with respect to the speed of light is 6t = 0.8 ± 0.7 sta t 

295 ± 2.9 sys ns, consistent with zero. The result obtained in 

296 Oct. -Nov. 201 1 is also consistent with zero, but it is sig- 

297 nificantly less precise. Taking into account the fact that 

298 some of the errors affecting the two measurements are 

299 correlated (the CERN side is the same and some Borex- 

300 ino detector systematics are also the same) we have de- 

301 cided not to average the two numbers and keep May 

302 2012 data only for the final result. 

303 The statistical error comes from the fit to the distribu- 

304 tion shown in Fig. [5] The systematic error depends on 

305 many sources that are listed and quantified in the Table 

306 below. 



307 Conclusions 

308 We have measured the speed of CNGS muon neutri- 

309 nos with the Borexino detector. 

Setting D = \x hx - £ C fiwl=730472.082 + 0.038 m as 
the distance between the target at CERN and the center 
of the Borexino detector, the final result for the relative 
time-of-flight distance between <E>=17 GeV muon 
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Description 


Error (ns) 


Time-Link Calibration (GPS) 


1.1 


Borexino electronics delays 


0.5 


Delays at CERN 


2.2 


Light propagation in BX detector 


1.0 


Electronics resolution 


0.5 


Event selection stability 


1.0 


Geodesy measurement 


0.1 


Total systematic error 


2.9 



Table 1: Table of systematic errors contributing to the May 2012 mea- 
surement. 



neutrinos with respect to the speed of light, 

D 

St - t BX - tcERN 

C 

310 is Si = 0.8 ± 0.7 + 2.9 ns, consistent with zero. From 

311 this number we obtain |v-c|/c < 2.1 10~ 6 90% C.L.. 
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